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Traditional calculations in perturbative quantum chromodynamics (pQCD) are based on an order- 
by-order expansion in the strong coupling as- Observables that are calculable in this way are 
known as “safe”. Recently, a class of unsafe observables was discovered that do not have a valid as 
expansion but are nevertheless calculable in pQCD using all-orders resummation. These observables 
are called “Sudakov safe” since singularities at each order are regulated by an all-orders Sudakov 
form factor. In this paper, we give a concrete definition of Sudakov safety based on conditional 
probability distributions, and we study a one-parameter family of momentum sharing observables 
that interpolate between the safe and unsafe regimes. The boundary between these regimes is 
particularly interesting, as the resulting distribution can be understood as the ultraviolet fixed 
point of a generalized fragmentation function, yielding a leading behavior that is independent of as- 


Infrared and collinear (IRC) safety has long been a 
guiding principle for determining which observables are 
calculable using perturbative quantum chromodynamics 
(pQCD) [Ilia. IRC safe observables are insensitive to 
arbitrarily soft gluon emissions and arbitrarily collinear 
parton splittings. This property ensures that perturba¬ 
tive singularities cancel between real and virtual emis¬ 
sions, leading to finite cross sections order-by-order in 
the strong coupling At the Large Hadron Collider 
(LHC), IRC safe jet algorithms like anti-fer play a 
key role in almost every analysis, and many jet-related 
cross sections have been calculated to next-to-leading 
and even next-to-next-to-leading order [M]. Of course, 
there are observables relevant for collider physics that are 
not IRC safe, though one can often use non-perturbative 
objects—like parton distribution functions, fragmenta¬ 
tion functions (FFs), and their generalizations [MU] — to 
absorb singularities and restore calculational control. 

In this paper, we show how to extend the calculational 
power of pQCD into the IRC unsafe regime using purely 
perturbative techniques. We study a class of unsafe ob¬ 
servables that are not defined at any fixed order in as, 
yet nevertheless have finite cross sections when all-orders 
effects are included. These observables are known in the 
literature as “Sudakov safe” m, since a perturbative 
Sudakov form factor m naturally (and exponentially) 
regulates real and virtual infrared (IR) divergences. To 
date, however, the study of Sudakov safe observables has 
been limited to specific examples. Here, we achieve a 
deeper understanding of these observables by providing 
a concrete definition of Sudakov safety based on condi¬ 
tional probabilities. The techniques in this paper apply 
to any perturbative quantum field theory, but we focus 
on pQCD to highlight an example of direct relevance to 
jet physics at the LHC. 

Because Sudakov safe observables are not defined at 
any fixed perturbative order, they in general have non- 
analytic dependence on as- Examples in the literature 
include observables with an apparent expansion in 
m and observables which are independent of 8.t suffi¬ 


ciently high energies [Ill [15]. As a case study, we consider 
a one-parameter family of momentum sharing observ¬ 
ables Zg based on “soft drop declustering” na, which al¬ 
ready appears in many jet substructure studies, e.g. nsi- 
ITH] . This family not only interpolates between the above 
two Sudakov-safe behaviors but also includes an IRC-safe 
regime. We explain how the boundary between the safe 
and unsafe regimes can be understood using the more 
familiar language of (generalized) FFs; the renormaliza¬ 
tion group (RG) evolution of the FF has an ultraviolet 
(UV) fixed point, suggesting an extended definition of 
IRC safety. 

To begin our general discussion of Sudakov safety, con¬ 
sider an IRC unsafe observable u and a companion IRC 
safe observable s. The observable s is chosen such that 
its measured value regulates all singularities of u. That 
is, even though the probability of measuring u, 


p{u) 


1 do¬ 
er drt ’ 


( 1 ) 


is ill-defined at any fixed perturbative order, the proba¬ 
bility of measuring u given s, p{u\s), is finite at all per¬ 
turbative orders, except possibly at isolated values of s; 
e.g., s = 0. Given this companion observable s, we want 
to know whether p{u) can be calculated from pQCD. 

Because s is IRC safe, p{s) is well-defined at all pertur¬ 
bative orders (although resummation may be required to 
regulate isolated singularities, see below). This allows us 
to define the joint probability distribution 


p{s,u) = p{s)p{u\s), (2) 


which is also finite at all perturbative orders, except pos¬ 
sibly at isolated values of s. To calculate p{u), we can 
simply marginalize over s: 


p{u) 


J dsp{s) p{u\s). 


( 3 ) 


Ifp(s) regulates all (isolated) singularities ofp{u\s), thus 
ensuring that the above integral is finite, then we deem 
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u to be Sudakov safe. In the case that one IRC safe 
observable is insufficient to regulate all singularities in 
u, we can measure a vector of IRC safe observables 
s = This gives a more general definition 

of Sudakov safety: 

p{u) = j d"'sp(s)p(u|s). (4) 

All previous examples of Sudakov safety fall in the cat¬ 
egory of ^ above where only a single IRC safe compan- 
sion s was required. In [14], the energy loss distribu¬ 
tion from soft drop grooming was dehned precisely as in 
(§, where u was the factional energy loss A^; and s was 
the groomed jet radius Vg (see below). In [T2] , ratio ob¬ 
servables r = a/h were originally dehned in terms of a 
double-differential cross section [191 ED] as 

p(r) = J dadbp{a, b) 6 (r - , (5) 

where a and b are IRC safe but r is not, because there are 
singularities at & = 0 at every hnite perturbative order, 
leading to a divide-by-zero issue for r. Integrating over 
a and using the dehnition of conditional probability ([^, 
we can write ([^ as 

p{r) = J dbp{b)p{r\b), (6) 

and r is Sudakov safe because p{b) has an all-orders Su¬ 
dakov form factor that renders p{r) hnite. 

It should be stressed that the dehnition of Sudakov 
safety in Q is not vacuous and it does not save all IRC 
unsafe observables. As a counterexample, consider par¬ 
ticle multiplicity; because perturbation theory allows an 
arbitrary number of soft or collinear emissions, one would 
need to measure an inhnite number of IRC safe observ¬ 
ables to regulate all singularities to all orders. Also, it 
should be stressed that just because an observable is Su¬ 
dakov safe, that does not imply that non-perturbative 
aspects of QCD are automatically suppressed. While a 
detailed discussion is beyond the scope of this paper, both 
[min] include an estimate of non-perturbative effects, 
which are analogous to power corrections and underly¬ 
ing event corrections familiar from the IRC safe case. In 
some cases, these corrections are known to scale away as 
a (fractional) inverse power of the collision energy. 

Crucially, one needs some kind of all-orders informa¬ 
tion to obtain hnite distributions for p{u). If a hxed- 
order expansion of p(s) and p{u\s) were sufficient, then 
p{u) would have a series expansion in a^, contradicting 
the assumption that u is IRC unsafe. In this paper, we 
use logarithmic resummation to capture all-orders infor¬ 
mation about p{s), which regulates isolated singularities 
at s = 0 to ensure the integral in (|^ is hnite. In all cases 
we have encountered, a hnite p{u\s) with a resummed 
p{s) is sufficient to calculate p(u), though this may not 
be the case generally. 


Unlike IRC safe distributions which have a unique as 
expansion, the formal perturbative accuracy of a Sudakov 
safe distribution is potentially ambiguous. First, there 
are different choices for s that can regulate the singular¬ 
ities in u. This is analogous to the choice of evolution 
variables in a parton shower, as each choice gives a hnite 
(albeit different) answer at a given perturbative accuracy. 
Second, the probability distributions p{s) and p(ri|s) can 
be calculated to different formal accuracies. Below we 
use leading logarithmic resummation for p{s), but only 
work to lowest order in as for p(u\s). Thus, when dis¬ 
cussing the accuracy of p{u), one must specify the choice 
of s and the accuracy of p{s) and p(u|s) separately. We 
stress, however, that the accuracy of both objects is sys¬ 
tematically improvable. 

We now study an instructive example that demon¬ 
strates the complementarity of Sudakov safety and IRC 
safety. This example is based on soft drop declustering 
[14] . which we briehy review. Consider a jet clustered 
with the Cambridge-Aachen (C/A) algorithm [51] 
with jet radius Rq. One can decluster through the jet’s 
branching history, grooming away the softer branch until 
one finds a branch that satisfies the condition 


min(pTi,PT2) ^ fR i2\ 

> ^cut Ip)’ ( ‘ ) 

PTI + PT2 \ Rq j 

where 1 and 2 denote the branches at that step in 
the clustering, pxi are the corresponding transverse mo¬ 
menta, and i?i 2 is their rapidity-azimuth separation. The 
kinematics of this branch defines the groomed jet radius 
Tg and the groomed momentum sharing Zg, 


Ri 2 mm{pTi,PT2) 

^3 = IT’ ^3 = „ , „ —; ( 8 ) 

Tg is IRC safe and its distribution was studied in [14] . 

Our observable of interest is Zg, and the angular ex¬ 
ponent /3 determines whether or not Zg is IRC safe. For 
13 < 0, Zg is IRC safe, because Zg > ^cut for any branch 
that passes Q; if this condition is never satished, the 
jet is simply removed from the analysis. For /3 > 0, 
Zg is IRC unsafe, since measuring Zg does not regulate 
collinear singularities. The boundary case /3 = 0 corre¬ 
sponds to the (modihed) mass drop tagger [TMT5] which 
also has collinear divergences, but we will show that it 
actually satisfies an extended version of IRC safety. 

In our calculations, we work to lowest non-trivial order 
to illustrate the physics, though we provide supplemen¬ 
tal materials for the interested reader that include higher- 
order (and non-perturbative) effects. We take the param¬ 
eter Zcut to be small, but large enough that log Zcut terms 
need not be resummed, with a benchmark of Zcut — 0.1. 

We now use the strategy in ^ to calculate the mo¬ 
mentum sharing Zg for all values of /3, using the groomed 
radius Tg to regulate collinear singularities: 

Pizg) = J ^'^9Pi'^g)p(zg\rg). (9) 
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FIG. 1. Distributions of Zg for various /3 values, obtained from 
Q at fixed Os = 0.1 and Zcut = 0.1. 


We use all-orders resummation to determine p{rg) and 
regulate the isolated = 0 singularity. This has been 
carried out to next-to-leading-logarithmic accuracy in 
M- Here, it is sufficient to consider the fixed-coupling 
limit: 


2asC^ d9 

TT X, ^ 

( 10 ) 

where Ci is the color factor of the jet, Pi{z) is the appro¬ 
priate splitting function (summed over final states), and 
the phase space cut is 


Ph) = g-exp 


[ dzP^{z)<^cnt 
Jo 


0cut = 0(1/2 - z)Q {z - ZcutO^) 

+ Q{z- 1/2) 0 ((1 - z) - . (11) 

The exponential part of ( |Io| ) is the rg Sudakov form fac¬ 
tor, where 0cut defines the no-emission criteria. To calcu¬ 
late p{zg\rg), note that Zg is defined by a single emission 
in the jet. For small Rq, the lowest-order matrix element 
is well-approximated by a 1 —>■ 2 splitting function: 


Pizg\rg) = ■' - 0(^3 - Zcutr^g ), ( 12 ) 

•^cut' g 

where 0 < Zg < 1/2 and we have introduced the notation 
P,(z) =P,(z)+P,(l-z). (13) 


In the double-logarithmic limit, we simply have Pi{z) = 
z^ allowing an explicit evaluation of 


P{Zg) = 


OisCi 


exp 


7r/3 


log 


(erf y 


a,Ci 

7r/3 




— erf 


P^iZg) (14) 


where 


/3 > 0 : oi = 0, 02 = min [2zcut, 2zg], (15) 

P <0 : ai = 2zg, 02 = 2zcut- (16) 



Safety 

Divergences 

Expansion 

d < 0 

IRC 

None 

„ Tl 

Qlg 

/3 = 0 

IRC via FF 

Collinear Only 

_ n — 1 

CXs 

d > 0 

Sudakov 

Collinear & Soft-Coil. 

nil 

as 


TABLE I. As d is adjusted, p(zg) interpolates between IRC- 
safe and two Sudakov-safe behaviors, related to the diver¬ 
gences in Zg. Here, n > 1 ranges over positive integers. 


Because ( [m] ) is finite, we see that Zg is at least Sudakov 
safe for all (3. Distributions of Zg calculated with at 
fixed as are shown in Fig. 

By expanding p{zg) in small Os, we can better under¬ 
stand the difference between IRC-safe and Sudakov-safe 
behavior. For /3 < 0, Zg is IRC safe, so Zg should have a 
well-defined expansion in as- To the accuracy calculated, 
([^ is fully valid to 0{as) in the collinear limit, and the 
expansion of ([^ yields the expected IRC safe result: 

/? < 0 : p(zg) = log —“ ^cut) 

'^\p\ ^cut 

+ Oia^s). (17) 

For j3 > 0, Zg is only Sudakov safe and its distribution 
should not have a valid Taylor series in Og. Indeed, for 
/3 > 0, the distribution has the expansion 


/3 > 0 : p{zg) = J'^^Pi{zg) + Oias), (18) 


and the presence of implies non-analytic dependence 
on Qfg. To 0{^/a^), the only phase space constraint is 
0 < Zg < 1/2, and the kink visible in Fig. 0 Zg — .2:cut 
first appears at 0{as). Finally, for the boundary case 
/I = 0, p{zg\rg) is independent of rg (in the fixed-coupling 


approximation), and (14) is independent of as'. 


P = 0: p{zg) = 


P^{Zg) 




^ -0(Zg - Zcut) ■ (19) 


We will later show that the /3 = 0 case does have a 
valid perturbative expansion in Og, despite being Og- 
independent at lowest order. The behavior of Zg for dif¬ 
ferent /? values is summarized in Table |Tj 

The /3 = 0 distribution of Zg is fascinating (and sim¬ 
pler than previous Og-independent examples [Til [IS]). 
Because Zg only has collinear divergences, we can un¬ 
derstand p{zg) in a different and illuminating way using 
FFs. As is well known, FFs absorb collinear divergences 
in final-state parton evolution, and we can introduce a 
generalized FF, F{zg), to play the same role for Zg. In 
the standard case, FFs are non-perturbative objects with 
perturbative RG evolution. In the Zg case, F{zg) is still 
a non-perturbative object, but it has a perturbative UV 
fixed point, becoming independent of IR boundary con¬ 
ditions at sufficiently high energies. 
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At Born level, the jet has a single parton, so Zg is 
undefined. We can, however, define F[zg) to be the one- 
prong Zg distribution, such that F{zg) acts like a non¬ 
trivial measurement function that is independent of the 
kinematics. Working to 0{as) in the collinear limit. 


P{Zg) = F{Zg) + 


UsCj ^ 

TT In 6 


X (Pi{Zg)Q{Zg-Zcut)-F{Zg) ^ (Iz Pi{z) 

\ J 2cut y 


+ C’(ag). 


( 20 ) 



There are two terms at 0{as). The first term accounts 
for the resolved case where the jet is composed of two 
prongs from a 1 —> 2 splitting. The second term corre¬ 
sponds to additional one-prong configurations (with the 
same F{zg) measurement function as the Born case), aris¬ 
ing either because the other prong has been removed by 
soft drop grooming or from one-prong virtual corrections. 

For a general F{zg), (20) is manifestly collinearly di¬ 
vergent because of the 6 integral, and F{zg) must be 
renormalized. But there is a unique choice of F{zg) for 
which collinear divergences are absent (at this order), 
without requiring renormalization: 


FMz,) = e(^g - ^cut). (21) 

J dzPi{z) 


Plugging this into ( [20| ), the 0{as) term vanishes, and we 
recover precisely the distribution in (19). 


In this way, Zg at (3 = 0 exhibits an extended version of 
IRC safety, where a non-trivial (and finite) measurement 
function is introduced in a region of phase space where 
the measurement would be otherwise undefined. Similar 
measurement functions appeared (without discussion) in 
the early days of jet physics [531 [21], where symmetries 
determined their form. Here, we used the cancellation 
of collinear divergences order-by-order in as to find an 


appropriate F{zg). We can also extend (20) beyond the 


collinear limit by considering full real and virtual matrix 
elements, leading to finite 0{as) corrections to p{zg). 

As alluded to above, Fuv(-Zg) also has the interpre¬ 
tation of being a UV fixed point from RG evolution. 
The collinear divergence of ( [20| can be absorbed into a 
renormalized FF, /x), at the price of introduc¬ 

ing explicit dependence on the MS renormalization scale 


/i. Requiring (20) to be independent of /x through 0{as) 
results in the following RG equation for F^^™\zg] p): 


a/i TT 


(— 

X P^{Zg)Q{Zg - Zeut) - F^'^-^\Zg-p) / , 

\ -'Zcut 


( 22 ) 


dzPi{z) . 


As p goes to -boo, the IR boundary condition is sup¬ 
pressed and F^’^‘^''\zg] fi) asymptotes to Fuy{zg). 


FIG. 2. Distributions of Zg for (3 = 0 and z^ut = 0.1 at the 
13 TeV LHC, as simulated by HerwigH —|- 2.6.3. The pr of 
the jets ranges from 50 GeV to 2 TeV, and the asymptotic 
distribution for quark jets, F((y in (211, is solid black. 


This UV asymptotic behavior can be tested using par- 
ton shower Monte Garlo generators. In Fig. [^ we show 
the Zg distribution for /3 = 0 for Herwig-|--|- 2.6.3 [53] 
at the 13 TeV LHG, using FastJet 3.1 [5^ and the 
RecursiveTools contrib [27]. As shown in the supple¬ 
ment, other generators give similar results. As the jet pt 
increases, p{zg) asymptotes to the form in (21) (which 
happens to be nearly identical for quark and gluon jets). 
This is due both to the RG flow in (22), which suppresses 
non-perturbative corrections, and the decrease of as with 
energy, which suppresses 0{as) corrections to p{zg). 


In this paper, we gave a concrete definition of Su- 
dakov safety, which extends the reach of pQGD beyond 
the traditional domain of IRG safe observables. Even 
at lowest perturbative order, the Zg example highlights 
the different analytic structures possible in the Sudakov 
safe regime, and the FF approach to the IRG safe/unsafe 
boundary yields new insights into the structure of per¬ 
turbative singularities. In addition to being an interest¬ 
ing conceptual result in perturbative field theory, ([^ of¬ 
fers a concrete prescription for how to leverage the grow¬ 
ing catalog of high-accuracy pQGD calculations (both 
fixed-order and resummed) to make predictions in the 
IRG unsafe regime. This can be done without have to 
rely (solely) on non-perturbative modeling, enhancing 
the prospects for precision jet physics in the LHG era. 
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SUPPLEMENTAL MATERIAL 

The supplemental material contains five additional cal¬ 
culations and analyses to highlight the behavior of Su¬ 
dakov safe observables, especially at higher orders. 


F. Krauss, M. Schon- 
JHEP 0902, 007 (2009) 


Calculating Ratio Observables 


When calculating the ratio observables in [T^] , there is 
a subtle difference between ([^ and § with respect to 
resummation. In [unn] which followed ([^, a Sudakov 
form factor for the joint probability distribution p{a, b) 
was found, which resums (some) logs of a, b, and a/b. In 
contrast, (|^ suggests first resuming logs of b in p{b), and 
then (optionally) resuming logs of r in p(r\b). These give 
slightly different expressions, as we show here. 

As in mus], the observables a and b are recoil-free 
angularities measured on a jet, which are defined as 




1 

Ptj 


iGJ 



(S-1) 


for angular exponent a > 0 and jet radius Rq. is the 
distance in the pseudorapidity-azimuth plane between 
particle i and the recoil-free broadening axis b [28] . 

For two angularities Cq, and ep, with a > j5, the joint 
cumulative probability distribution at leading logarith¬ 
mic accuracy with fixed coupling is m 


= exp 


f log^ 6/3 
TT * y P 


a-P j 


(S-2) 


where Ci is the color factor. The double differential cross 
section/joint probability distribution is then 


P{ea,ep) 


1 dV 

a Aca Aep 


92 


dCa dCjS 


Xj(eQ,, e^). 


(S-3) 


To determine the distribution for the ratio r = e^/e^, 
we can then insert this into ([^. The full expression is 
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given in and the lowest order terms in the small as 
expansion are 


p{r) 


1 da 
a dr 



v^l 

a — P r 


de^dep ^ J S (r- 
dsct dey3 y 

_ ^as Ci(a - 2/3) logr 

TT (a — PY r 



+ 0[aT). 


This distribution can be considered the leading logarith¬ 
mic distribution for r because it was calculated from the 
double differential cross section calculated to leading log¬ 
arithmic accuracy in a, 6, and a/h. 

Alternatively, we can calculate p{r) as in Q, where 
the conditional probability, calculated to fixed order, is 
integrated against the resummed distribution for the de¬ 
nominator. To leading logarithmic accuracy with fixed 
coupling, the cumulative distribution for ep is 


T.{ep) = exp 


TT P 



(S-5) 


The probability distribution of ep is thus 


Calculating the conditional probability using the most 
singular terms in the splitting function, we find 


p{r\ep) 



1/r 

P log 1/e^ 


0(1-r)0 



(S-7) 


to lowest order, which is indeed a normalized conditional 
probability distribution. Using the method of yields 
the probability distribution for the ratio r: 


where /(s) is some function, and E(s) = is the 

resummed cumulative distribution for s. Depending on 
whether s has both soft and collinear singularities or just 
collinear ones, /(s) is expected to take different forms: 

OO 

fsc{s) = ^ Cfc s + ..., (S-10) 

A .-1 

OO 

fc{s) = '^dka';\og''s + . ■., (S-11) 

where the ellipses (...) stand for terms with additional 
as suppression. We say that fsc{s) has “double log” be¬ 
havior (since the lowest term is as log^ s), while /c(s) is 
“single log” (for a^logs). 

The functional form of p{u\s) is not needed to derive 
our heuristic, though a few facts about p{u\s) are im¬ 
portant. First, since s regulates the divergences in u, 
p{u\s) must have a valid Taylor expansion in as- Second, 
because s itself has singularities, p(m|s) will necessarily 
have dependence on log s related to the structure of p(s). 
Third, because the conditional probability distribution 
is normalized as f dup(s|u) = 1, p(u\s) has to start at 
0{a^). In particular, at lowest order 

Psc(u|s) = ^f^+0(a,), (S-12) 

logs 

Pc{u\s) = h{u,s) + 0{as), (S-13) 

where we have pulled out an extra log s factor in the 
soft/collinear case, such that g{u, s) and h{u, s) only have 
power-suppressed dependence on s. 

From these generic forms for p(s) and p(u\s), we can 
determine p{u) using ([^, 


p{r) = J dep p{ep) p{r\ep) (S-8) 

.^/PCll , CiP logr ^ ^^^ 3 / 2 ^ 


= va, 


a — P r TT {a — PY 


The term at agrees with (S-4), but the term 

at 0{as), and higher terms, generically do not. This 
emphasizes that to determine the formal accuracy of a 
Sudakov safe observable requires specifying the accuracy 
of all components in its calculation. 


Heuristic for the Os Expansion 

In Table [ij we drew a link between the singularities 
present in Zg and the expected as expansion for p{zg). 
Here, we give a heuristic way to understand this behavior. 

For a generic resummed observable s G (0,1), the prob¬ 
ability distribution p(s) can be written as 


p{u) = j ds f'{s) p{u\s) (S-I4) 

= J dx p{u\s{x)), (S-I5) 


where we have introduced the change of variables 

x = f{s), dx = f{x)ds, (S-I6) 


and we expect x G (—oo,0). With this change of 
variables, the as dependence of p{u) resides entirely in 
p(u\s{x)). In general, there is no closed form for s{x), but 
we can determine it order by order in as by inverting the 
series in (S-10) and (S-111. Depending on the divergence 
structure of s, there are different forms: 


sc: as log s(a:) 

c : as logs(a:) 


y^ + 0(as), 

^TC^(as). 


(S-17) 

(S-18) 


p{s) 


Ae/(«) - 
ds “ 


f{s) e 


fis) 


(S-9) 


In the soft/collinear case, the y^as factors in logs mean 
that p(M|s(a;)) will have a expansion when expressed 
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as a function of x. By contrast, in the collinear only case, 
p(M|s(a:)) will still have an ordinary expansion. 

The last ingredient is the starting order of the expan¬ 
sion. Plugging in the lowest order expressions for p(u|s), 
we can evaluate p(u) using (S-15), up to power correc¬ 
tions: 


Psciu) = V-asCiTTg{u, 0) -I- C>(as), (S-19) 

Pc{u) = h{u, 0 ) + 0 {as). (S-20) 

This confirms the expected expansions in Table |lj Note 
that in the collinear only case, there is no leading depen¬ 
dence on di, though di will show up at 0{as). 


The Zg Distribution at Higher Accuracy 


Via Q, we calculated p{zg) to the lowest non-trivial 
order. Here, we discuss how to improve the accuracy of 
this calculation through running coupling effects. The 
following discussion is valid for any f3. 

We start with the groomed radius distribution pirg), 
for which the all-orders resummation was derived in El- 
Including running coupling effects in (10), we obtain 




P{fg) = 

where 



II 

2a 

de f 


Lg 0 Jo 


_2a 

de r 


TT 

Irg T ^ 


dr„^ 


f(rg) 


(S-21) 


dzPt(z) as {zBptRo) © cut 

^ 1/2 _ 

= 1^1 dzPi{z)as{z0pTRQ), 




(S-22) 


z = min( 2 :, 1 — z), and ©cut is defined in 0- When 
summed over final states, the quark and gluon splitting 
functions are 


P.iz) = 

PM = 


l + il-zf 
2z 

1-z z{l-z) 




(S-23) 


(S-24) 


Exploiting z (1 — z) symmetry, we have written the 
gluon splitting function in such a way that it exhibits a 
singularity only when z —>■ 0. As noted in [2], (S-211 is 
accurate to single-logarithmic accuracy, provided that as 
is evaluated in the CMW scheme [35]. This is because 
Xg is set by just one splitting, with no multiple-emissions 


contribution. The expression in (S-211 only differs from 


the corresponding one in El because of the more sophis¬ 
ticated ©cut treatment, which takes into account hnite z 
corrections. 



FIG. S3. Resummed distributions of Zg for Zcut = 0.1 and 
three different values of /3 = 0.5, 0, —0.5. We take pr = 
2 TeV and Ro = 0.5. Solid curves represent resummation 
with running coupling, while dashed ones are evaluated at 
fixed coupling as{pTRo) = 0.087. 


Next, we address the conditional probability p{zg\rg). 
There are various strategies to compute this quantity. 
For example, we could use exact hxed-order matrix ele¬ 
ments or we could embark on a systematic all-orders cal¬ 
culation. Here, we will show an intermediate approach, 
working in the collinear limit for the matrix element, but 
including a tower of all-orders contributions originating 
from as running. Specihcally, we write 


PiMg) 


P*i^gPg) 

P*irg) 


©(zg 



(S-25) 


where 


P*{zg,rg) 


2a /-i 

- / TT / dz Pi[z)<dcutO{z - Zg) 

Jo 0 Jo 

X d{0 - Tg) as {zOptRo) 

2C _ 

— -asizgXgPTRo) Pi{zg) (S-26) 


has as evaluated at the proper scale of the emission, and 



dZgP*{Zg,rg) 




dzgPi{zg)as 


(ZgrgPTRo) 


(S-27) 


ensures that p{zg\rg) is properly normalized. 

With these running coupling improvements, p(zg) can 
be computed using the Sudakov safe definition in (|^. To 
regularize the Landau pole in as, we freeze its running 
at the non-perturbative scale = 1 GeV. Our results 
are shown in Fig. |S3| and compared to the fixed coupling 
case with asippRo) = 0.087. 
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Renormalization Group Evolution of the 
Fragmentation Function 


which then produces the homogeneous solution 

dzPi(z) 


As discussed in the paper, the collinear divergence of (^91 m) Fo{zg,fJ,o) 


(20 1 can be absorbed into a renormalized FF, at the price 






of introducing explicit scale dependence. In dimensional 
regularization in the MS scheme, the renormalized F(zg) 
to 0{as) is 




sC, 


(S-28) 


r.1/2 


X \ Pi{Zg)Q{Zg - Zcut) - F{Zg) / dzP^{z)] , 


where e is the dimensional regularization parameter, 
c = log( 47 re“'''®), and /r is the renormalization scale. Re¬ 
quiring the cross section to be independent of ^ through 
0{as) results in an RG equation for {zg-, fj,): 


7T 


(S-29) 


X (pi{zg)Qizg - Zcnt) - F'^'^^''\zg;fi) f dzPi(z) j , 
\ ^cut / 


up to corrections at 0{a^). Compared to (22 1 , here we 
have explicitly introduced the /i dependence in as- 
To solve for {zg; ^), we first find the solution to 

the homogeneous equation: 


^ j-,(ren)/ \ 




dz Pi{z) 


The solution is 

= ^ 0 ( 2 ^ 9 ; Mo) exp 


(S-30) 

(S-31) 


^ J Zcut 9 Un M 


where FQ{zg; ^q) is the boundary value defined at the 
infrared scale fj,Q. Using the definition of the /3-function, 


a/ \ 


(S-32) 


the integral over the scale /r' can be exchanged for an 
integral over as itself: 


rdM' . 

/ 

Vo ^ 


sip) da' 


'as(/zo) ) 

The one loop ,0-function is 

Picxs) = -alPo , 


a . 


(S-33) 


(S-34) 


(S-35) 

The exponent is found by integrating the quark and gluon 


splitting functions in (S-23) and (S-24): 
.1/2 


' Zcut 

. 1/2 


dz Pq{z) = log ( 


dzPg{z) = log 


1 -^cut 


-k 


uf _ n\ 


Zcnt J \dCA 12/ 

/ /. 4-' V X. 

22 :cut 1 1 


_ ^ f 1 _ ZLh 1 

dCAj 2 I Ca 


^cut f 1 _ 

3 I CaJ' 


(S-37) 


Now, we must find a particular solution to ( |S-29 ). The 
simplest approach is to assume that Fd^’^)(zg; fj,) is inde¬ 
pendent of /i, which requires: 

e(z, - Zeut). (S-38) 

iz., 


The full solution to ( |S-29 | is then the sum of the homo¬ 
geneous and particular solutions: 


= ^ 0 (^ 9 ; ^ 0 ) 


p \~9 

as{n) \ ^ 
a!s(Mo) 


pp dzPi(z) 


Fi{Zg) 


- 0(^3 -Zcut). (S-39) 


By asymptotic freedom of cis, as /r —>■ 00 , the homoge¬ 
neous solution is suppressed and any dependence on the 
infrared boundary condition Fo(zg;/io) vanishes. There¬ 
fore, the asymptotic distribution is 

Fd“)(z9; fj,^oo) = - y V ~ - -“ ^cut) • 

(S-40) 


This agrees with the Sudakov safe calculation in (19), and 


makes a definite prediction for the high energy behavior 
of this observable. 


We can also use (S-391 to estimate the scaling of non- 
perturbative corrections to p{zg). The boundary condi¬ 
tion Fo{zg; pq) is 0 ( 1 ) and fully non-perturbative, but in 
the small Zcut limit it scales away approximately as 


l-tas(Mo)/3olog — 
Mo 


(S-41) 
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Though not quite a power law suppression, it would be 
a power law in the /3o —t 0 limit 


■ log : 


■0(/3o), 


(S-42) 


with typical exponent ~ 0.1 (~ 0.2) in the quark (gluon) 
case. In general, because of the larger C^, we predict 
that gluon-initiated jets will saturate F\jy{zg) faster than 
quark-initiated jets. Beyond F[zg), there may be shape 
function corrections as with IRC safe observables. 

Finally, we have found evidence that the FF UV fixed 
point in (21) may be one-loop exact. We have verified 
this explicitly at 0(ag) in the strongly-ordered collinear 
limit but have not yet attempted a complete proof. We 
observe that the collinear singularity in Zg is resolved at 
any non-zero value of rg. We therefore expect that the 
only kinds of singularities present are those associated 
with a 1 —> 2 splitting, which are already included in 
(S-29). No additional collinear singularities are expected 
to appear in 1 —> 3 splittings, since either the third par- 
ton generates a finite value of Vg (in which case all sin¬ 
gularities are regulated), the third parton fails the soft 
drop condition (in which case the kinematics reverts to 
I —>■ 2), or the third parton is collinear with the other 
two (in which case the real I —3 collinear singularity 
should cancel against the virtual I —?► 2 one). 


Monte Carlo Analysis of Zg for /3 = 0 

When comparing our analytic calculation of p{zg) at 
/I = 0 to Monte Carlo generators in Fig. we only 
showed results for Herwig-|-+ 2.6.3 [2S]. Here, we in¬ 
clude two additional generators: Pythia 8.201 [30] and 
Sherpa 2.1.1 |3T|. Although not reported here, we 
checked that the features below are also observed in the 
ViNCiA 1.1.3 |32| antenna shower applied to high en¬ 


ergy electron-positron collisions. In all cases, we used 
FastJet 3.1 |2^ to reconstruct jets and the Recursive- 
Tools contrib m to implement soft drop grooming. 

In Fig. we show the Zg distribution at /3 = 0 for 
three different values of the energy cut, Zcut = 0.2, O.I, 
and 0.05. All samples are at the 13 TeV LHC, includ¬ 
ing hadronization effects as well as the default underly¬ 
ing event models. For simplicity we only show the Zg 
distributions for jets with large transverse momentum, 
Pt > 2 TeV, in order to reduce hadronization correc¬ 
tions. Jets at Zg = 0 correspond to situations where the 
soft drop groomer gives a one-prong configuration, and 
those events are not used for normalizing p{zg). 

All distributions are in decent agreement with the 
fixed-coupling analytic prediction in (21). The agreement 
is particularly good for larger values of Zcut but degrades 
as Zcut gets smaller, especially in the Pythia sample. 
We interpret the primary difference between our ana¬ 
lytic calculation and the Monte Carlo generators as aris¬ 
ing from running coupling effects (see Fig. S3) and from 
residual dependence on non-perturbative hadronization 
corrections. As Zcut decreases, the phase-space for emis¬ 
sions which build up the Sudakov form factor extends 
more and more into the soft region, eventually picking 
up contributions from non-perturbative emissions. This 
affects small values of Zg, which then distorts the overall 
p{zg) distribution after normalization. 

Though not shown, the small Zcut distortion is more 
pronounced for quark jets than for gluon jets. This is 
expected from the analysis of (S-4I), since the larger 
gluon color factor means that more of the distribution 
is described by wide angle perturbative emissions. An¬ 
other potential source of distortion is large logarithms 
of Zg and Zcut which have not been resummed. Some of 
these logarithms can be captured through running cou¬ 
pling effects, and preliminary investigations indicate that 
remaining logarithms may give a small contribution, ow¬ 
ing to the single-emission nature of Zg. 
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Groomed Momentum Fraction 


13 TeV LHC, Pt> 2 TeV 


Ro = 0.5, =0, zcu, 
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FIG. S4. Distributions of Zg for /I = 0 at the 13 TeV LHC as 
simulated from HerwigH— h, Pythia, and Sherpa Samples. 
From top to bottom, the Zcut values are 0.2, 0.1 and 0.05. 
Jets are selected with pt > 2 TeV and jet radius Rq = 0.5. 


















